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Abstract 



■^ , The low-temperature structure of ^Ceo (^=K, Rb) is an ordered array of 

^^ I polymerized Ceo chains, with magnetic properties that suggest a non-metallic 

._!. I ground state. We study the paramagnetic state of this phase using first- 

"^ I principles electronic-structure methods, and examine the magnetic fluctua- 

tions around this state using a model Hamiltonian. The electronic and mag- 
netic properties of even this polymerized phase remain strongly three dimen- 
. ,_, sional, and the magnetic fluctuations favor an unusual three-dimensional anti- 

H ferromagnetically ordered structure with a semi-metallic electronic spectrum. 
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The family of alkali-intercalated fullerides, AnC^o, has a remarkably rich phase diagram. 
For large alkalis (K, Rb, Cs, and their binary mixtures) a variety of stable phases have 
been identified, the electronic and magnetic properties of which are still poorly understood. 
Conventional band theory appears to explain the insulating n=0 and 6 terminal phases, and 
suggests a conventional metallic normal state for the n=3 phase, but it incorrectly predicts 
that the intermediate n=4 phase is also metallic. Various scenarios have been advanced to 
explain this discrepancy, primarily by appealing either to very strong intramolecular electron 
correlations ||l|J^ or to various possible spin- or charge-density-wave instabilities from the 
paramagnetic state. This latter point of view has received additional impetus from the recent 
discovery of a polymerized form of (n=l) RbCgo and KCgo [@]- Electron spin resonance data 
for this phase show a phase transition at 50 K, below which the spin susceptibility decreases 
by an order of magnitude, suggesting an insulating ground state. Chauvet et al. have 
proposed that the single alkali electron results in a half-filled conduction band, and that 
the polymer-like crystal structure then promotes the formation of a spin-density wave that 
doubles the unit cell along the polymer chain direction and opens a gap at the Fermi level 

i- 

In this Letter, we use the results from density-functional electronic-structure methods to 

construct a physically realistic model Hamiltonian to describe the electronic states in the low- 
temperature phase of ACqq. Remarkably, we find that the electronic and magnetic properties 
of even this polymerized phase remain strongly three dimensional. Indeed, by studying 
the fiuctuations around the paramagnetic reference state, we identify a broken-symmetry 
semi-metallic ground state with a three-dimensional magnetically ordered structure quite 
different from the quasi- ID structure previously proposed for this phase. We note, however, 
that the quasi-lD electronic character hypothesized by Chauvet et al. might be achieved 
by dilation of the galleries between polymer chains, through the use of larger intercalants. 
The transition from the three-dimensional magnetic structure discussed here to a quasi- ID 
structure is quite interesting and is beginning to be explored. 

ACqo {A=K,IVo) undergoes a first-order structural phase transition around 350 K, from a 



high-temperature fee phase to a low-temperature orthorhombic phase 0. A striking feature 
of the orthorhombic phase is the unusually short nearest-neighbor fuUerene distance of 9.12 A 
(for both K and Rb intercalants) [Q]. This intermolecular spacing implies a separation 
between interatomic carbon atoms on neighboring (icosahedral) fuUerenes as small as 2.0 A, 
leading to the suggestion that adjacent molecules are actually covalently bonded to form 
polymer-like chains (C6o)n aligned along one crystallographic axis |^. This picture has been 
confirmed by Rietveld refinement of x-ray powder diffraction data [Q . 

To study the electronic states in such a structure, we begin by performing electronic- 
structure calculations for the low-temperature orthorhombic phase of crystalline RbCeo, 
using the local-density approximation (LDA) to density-functional theory. For the details 
regarding the computational methods, see Ref. and references therein. We have made 
one structural simplification to the Rietveld model: by rotating every chain through 45° in 
the same sense, the Bravais lattice is changed from simple orthorhombic to body-centered 
orthorhombic (bco), and the number of formula units per unit cell is reduced from two to one. 
The resulting electronic conduction-band structure is shown in Fig. 1. Quite surprisingly, 
we find that in the occupied part of the spectrum, the dispersion perpendicular to the 
chains (F-M ) is quite comparable to the dispersion along the chains (T-H) . In other words, 
the amplitude for electron hopping is nearly isotropic, with the resulting Fermi surface 
exhibiting very strong three-dimensional character. Two other features of Fig. 1 deserve 
comment. First, the low symmetry of this structure admits only one-dimensional irreducible 
representations, so that the three-fold degeneracy of the parent ti„ state (of the cubic phase) 
is lifted. The magnitude of this symmetry breaking is large on the scale of the conduction 
band width, of order half a volt. Second, the lowest band, which is roughly half filled, shows 
positive dispersion at F, so that the occupied sector of the Brillouin zone encloses the zone 
center. We find that this lower band transforms like z, that is, these wave functions have 
odd parity under refiection through the mirror plane normal to the chains. 

Qualitatively, the comparable longitudinal and transverse dispersion is a consequence 
of two distinct structural features: (1) In order to covalently bond two adjacent Ceo 



molecules, the valence orbitals on intermolecular nearest-neighbor carbon atoms rehybridize 
from threefold-coordinated (radial tt) to fourfold-coordinated (approximately sp^) orbitals. 
Consequently, the 7r-like conduction band in this polymerized phase contains only a very 
small admixture of orbitals on intermolecular nearest-neighbor atoms, so that dispersion 
in this direction is actually governed by the interaction of more distant neighbors |^). (2) 
Even with relatively weak interchain electron hopping, dispersion across the chains also ben- 
efits from high molecular coordination in the transverse directions. With a larger interchain 
separation, this benefit would be largely lost and intrachain hopping would predominate, 
leading to a regime whose consequences are briefly discussed below. 

To study the stability of this paramagnetic reference state, we now focus our attention 
on the dynamics within the doped conduction band. We consider an effective Hamiltonian 
of the form 



H = Yl (4t,Tif,,juCju + h.C.) +UY, ^iMT^iMi' (1) 



J 

where c] creates an electron on site i with orbital polarization jj = x,y, z, and the Tj^ jv are a 
set of 3x3 matrices giving the amplitudes for an electron hopping between different orbitals 
on different sites. The electron-electron interaction terms describe a residual short-range 
repulsive potential which, for simplicity, we take to be diagonal in both the site and orbital 
indices. For the kinetic-energy term, we retain terms coupling each molecule to its first 12 
nearest neighbors; this closes the first coordination shell in the parent fee structure (the 
orthorhombic distortion breaks this into five inequivalent shells) . We also include a diagonal 
on-site term to correctly account for the orthorhombic crystal-field splitting of the three-fold 
tiu molecular state. The matrix elements Tj^^v are computed by direct Fourier inversion 
of the LDA eigenvalue spectrum; this generates the optimal tight-binding representation 
of the LDA bands, and closely reproduces the first-principles single-particle spectrum. (A 
detailed description of this parametrization procedure is given in Ref. |^.) The resulting 
tight-binding dispersion is shown in Fig. 1 as dotted curves; the quality of the fit to the 
LDA spectrum is clearly excellent. 



We now consider a mean-field decoupling of the interaction term in Eq. (|l]), and study 
the resulting (inverse) generalized density-wave susceptibility 

= -u {6,,6,j,,^^, + f/n°^.,^,,,) , (2) 

where Q is the grand potential, and 

and a labels an eigenfunction of the single-particle piece of Eq. (|l]) with eigenvalue E^ and 
occupation number /q. This interaction couples fluctuations in orbital fi on site i with those 
in orbital z/ on site j. Diagonalization of the susceptibility matrix then identifies both the 
wave-vector and orbital character of the dominant spin- and charge-density fluctuations in 
the model. We find that for a short-range repulsive potential, U > 0, the spin fluctuations 
within the 2;-polarized orbital are dominant. These fluctuations can be studied further by 
considering the matrix 

K'^Hq) = - E K k\P^a,{-q)\n, k + q) 

mnk 

x{n,k + q\aMP.\m, k) t"'" ~ ^"'^+'^ , (4) 

where P^ projects the /i-orbital polarization. By examining the momentum dependence of 
the largest eigenvalue, Hz{q), of the matrix K^^l{q) we find that the dominant spin-density 
fluctuations are quite strongly peaked near the X point of the bco Brillouin zone. (Note 
that the T-X direction is perpendicular to the chain direction, and the X point is the center 
of the zone face.) To illustrate. Fig. 2(a) shows the static spin structure factor 

^''^'^ = i-(?/2K(g)' ^^^ 

calculated at T=0 and f/=250 meV. This interaction strength is less than the zero- 
temperature critical value, t/c=265 meV, so that here the paramagnetic phase is stable 



at low temperature. The data are plotted in the qx-Qz plane, and the peaks of the correla- 
tion function are centered at the X point and its translates. Fig. 2(b) shows the same data 
as a linear plot along the T-X direction, calculated at three temperatures approaching the 
mean-field ordering temperature, Tc, for f/=280 meV. Within the mean-field theory, we find 
that a spin-density-wave transition at 100 K would require a repulsive potential U ~ 270 
meV, smaller than most current estimates of the effective intramolecular repulsion given by 
microscopic theories P,^. 

A condensation at the X point for T < Tc describes the ordering of a three-dimensional 
antiferromagnetic state with two sublattices, such that the spin polarization alternates be- 
tween the corner and body-centered sites of the bco structure. This spin configuration 
can thus be described as ferromagnetic within each chain and antiferromagnetic between 
nearest-neighbor chains. We note that the three-dimensional character of this magnetic 
state might have been anticipated in view of the substantial three-dimensional character 
of the single-particle spectrum displayed in Fig. 1. If one views the unpaired spins on the 
fullerene molecules as localized and antiferromagnetically coupled via Ji > to nearest 
neighbors along a cell diagonal, and antiferromagnetically via J2 > along nearest neigh- 
bors within the chains, then the transition to one- dimensional behavior requires J2/J1 > 2. 
This limit is apparently well outside the regime that actually describes these materials. In 
fact, we find that the effective intersite interaction mediated by the polarization of the con- 
duction electrons gives J2 < 0; that is, the nearest-neighbor intrachain interaction is actually 
ferromagnetic for this structure. 

Fig. 3 shows the spectral density obtained at T=0 in the equilibrium phase of this model; 
panel (a) is the density of states calculated for the unbroken paramagnetic phase, while panel 
(b) is for the equilibrium magnetically ordered phase. We find that for f/=300 meV the Fermi 
surface of the parent phase is partially gapped in the ground-state structure, and that a well 
developed pseudogap is clearly evident in the density of states. Thus the system remains 
semi-metallic at low temperatures. In the spectrum one can nevertheless clearly observe the 
signature of the ferromagnetic intrachain ordering, which essentially results in a rigid spin 



splitting of the lowest (z-polarized) conduction band, evident in the structure of the density 
of states around the Fermi level. The upper two unoccupied branches of the conduction 
manifold are relatively insensitive to this ordering transition. We find that the pseudogap is 
a robust feature of this ground state, and survives up to relatively large coupling strengths; 
a fully developed gap first opens in the single particle spectrum for ?7=470 meV. We suggest 
that the decrease in the spin susceptibility measured by Chauvet et al. may be associated 
with the opening of the pseudogap at T^. Further measurements at temperatures much lower 
than Tc will be required to distinguish between insulating and semi-metallic behavior in this 
phase. 

Our density-functional calculations indicate that the transverse interchain hopping am- 
plitudes in the Hamiltonian of Eq. (|l]) are very sensitive to the interchain separation. Indeed, 
we find that a 10% dilation of the interchain spacing leads to a regime where the electronic 
structure is dominated by one- dimensional dynamics along the chain direction. Here the spin 
fluctuations are again dominantly antiferromagnetic in character, but now with the ordering 
wave vector oriented along the chain. We emphasize that while this limit is apparently not 
reahzed in K- or Rb-doped samples, it might well be achieved by doping with substantially 
larger species. Experimentally, such a situation could conceivably be realized by alloying 
the alkali with "spacer molecules" such as NH3, as has been done for A3 Ceo [0- If so, the 
resulting magnetic phase diagram as a function of alloying composition is expected to be 
quite complex [|ll| . 



The results reported above have been obtained for a model structure with a common 
orientation on each of the bco Bravais-lattice sites, and therefore it is natural to ask how 
sensitive the results are to this assumption about the orientational order. While detailed 
calculations are continuing, we should note that for the Pmnn structure (determined by Ri- 
etveld refinement) with two inequivalent orientations per orthorhombic unit cell, our LDA 
calculations show that the low-energy spectrum retains its strong three-dimensional char- 
acter ||n[. Since this feature controls the physics discussed above, we expect the scenario 
developed above for the single-orientation structure still to apply. Also, we should note that 



the possibility of a quenched orientationally disordered phase for AC^o has not yet been 
excluded experimentally [Q. 

In summary, we have studied the paramagnetic phase of the polymerized fulleride ACqq 
and have examined magnetic fluctuations around this reference state. We find that repulsive 
interactions favor a two-sublattice three-dimensional magnetic structure for this system, 
which can be stabilized at relatively modest values of the repulsion strength. The three- 
dimensional character of this phase derives primarily from the rehybridization of the cubic 
tiu conduction band in this polymerized phase, and from the large coordination transverse 
to the polymer chain axis. Finally, we observe that a physically plausible modification of the 
crystal structure, possibly by intercalation with larger species, may be sufficient to stabilize 
the quasi-lD magnetic state which had been earlier hypothesized for this phase. 
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FIGURES 
FIG. 1. Theoretical LDA conduction band structure (solid curves) and tight-binding fit to 

the LDA spectrum (dotted curves) for polymerized RbiCeo, in the body-centered-orthorhombic 

structure described in the text. The labeling of symmetry points is taken from the more familiar 

body-centered tetragonal zone: H is the zone edge along the chain direction, and M is the zone 

edge in the perpendicular direction. The Fermi level is the energy zero. 

FIG. 2. (a) Surface plot of the static spin structure factor, SzziQ), for r=0 and U=250 meV, 
plotted in the Qx — Qz plane. The sharp peak at the X point of the bco zone corresponds to the 
three-dimensional antiferromagnetic ordering described in the text, (b) Szz{q) plotted along T-X, 
for U=280 meV and shown for three temperatures approaching the ordering value, Tc=130 K. 

FIG. 3. Spectral density at T=0 for (a) the paramagnetic reference state, and (b) the magnet- 
ically ordered ground state. The Fermi level is the energy zero. 
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Fig. 1 S.C. Erwin, G.V. Krishna, and E.J. Mele 




Fig. 2(a) S.C. Erwin, G.V. Krishna, and E.J. Mele 
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Fig. 2(b) S.C. Erwin, G.V. Krishna, and E.J. Mele 
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Fig. 3 S.C. Erwin, G.V. Krishna, and E.J. Mele 



